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Classical gravitational bremsstrahlung in particle collisions at transplanckian energies is studied in 
Mi X T"^. The radiation efficiency e = -Erad/£'initiai is computed in terms of the Schwarzschild radius 
'"s(\/s), the impact parameter b and the Lorentz factor 7cm and found to be e = Cdirs /b^'^^^'yl^^ , 
larger than previous estimates by many powers of 7cm 3> 1. The result is reliable for impact 
parameters in the overlap of < & < &c and 6 > \c, with b^ marking (for d 7^ 0) the loss of the 
notion of classical trajectories and Ac = h/mc the Compton length of the scattered particles. The 
condition on s and m for extreme radiation damping and (presumably) no black hole production is 
also derived. 



PACS numbers: 04.50.Cd, 04.50.Gh, 04.70.Dy 

Black hole (BH) production in LHC, predicted [l| 
in models with TeV-scale gravity and large extra di- 
mensions 0, m, Q about ten years ago, has been the 
subject of intense theoretical study and numerical sim- 
ulations (for a review see [H])- The prediction is based 
on the assumption that for impact parameters of the 
order of the horizon radius corresponding to the CM 
collision energy 2E = -^s 
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an event horizon should form due to the non-linear 
nature of gravity. The D = A + d dimensional grav- 
itational constant is Gd = h'^^^ /{M^+'^c'^-^), with 
the D-dimensional Planck mass. It is related 
to the four-dimensional Planck mass Mpi = M4 via 
Af2j = M^+'^V, V = {2tt RY, where R is the large 
compactification radius. 

This classical, essentially, picture of BH formation is 
justified for transplanckian energies s » G^}^'^'^^'^^ — 
Ml. Indeed 0, in this case the D-dimensional Planck 

length U = (SGd/c^)^/^'^"^^^ = hjKhc and the de 
Broglie length of the collision As = hc/^ satisfy 
the classicality condition Xb <C I* ^ ^5. Further- 
more, gravity is believed to be the dominant force in 
the transplanckian region. Thus, for BH masses large 
compared to M*, the use of classical Einstein theory 
is well justified. Moreover, it seems that formation 
of BHs in four dimensions is predicted by string the- 
ory Thus, in spite of the fact that there are is- 
sues which require further study 8], a consensus has 
been reached that the prediction of BHs in ultra-high 
energy collisions is robust and is summarized in the 
widely accepted four-stage process of formation and 
evaporation of BHs in colliders [1, 9i|, namely (i) for- 
mation of a closed trapped surface (CTS) in the col- 
lision of shock waves modeling the head-on particle 
collision, (it) the balding phase, during which the BH 
emits gravitational waves and relaxes to the Myers- 
Perry BH, (iii) Hawking evaporation and superradi- 
ance phase in which the experimental signatures are 
supposed to be produced, and (iv) the quantum grav- 



ity stage, where more fundamental theory like super- 
strings is important. This scenario was implemented 
in computer codes to simulate the BH events in 
LHC, where they are expected to be produced at a 
rate of several per second, and in ultra high energy 
cosmic rays. 

Here we focus on stage (i). Replacing the field of 
an ultrarelativistic particle by that of a black hole in 
the infinite momentum frame strictly speaking is only 
valid in the linearized level, while the associated non- 
linear phenomena for particles may be quite different 
from those in colliding waves. A typical non-linear 
effect for particles is gravitational bremsstrahlung, 
which is an important inelastic process in transplanck- 
ian collisions. The existing estimates of gravitational 
radiation either refer to phase (ii), or are based on 
the assumption of an already existing BH (e.g . ra- 
diation from particles falling into the BH [11|), on 
results of linearized theory relevant only to the case 
of non-gravitational scattering Q , on weakly relativis- 
tic numerical simulations [l^ or again on collisions of 
waves in 4d [l^l- However, a detailed study of gravi- 
tational bremsstrahlung in the transplanckian regime 
in the ADD scenario was missing. The purpose of this 
note is to present the results of such a study (a more 
detailed account will follow [l3|)- The emitted energy 
is found to be larger than earlier estimates (see e.g. 
@) by a dimension dependent power of the Lorentz 
factor of the collision. The condition is derived, for 
which radiation reaction [Tsj becomes essential and 
the picture of colliding waves, which models particle 
motion with constant speed, may not be adequate. 

The standard ADD model assumes that empty 
space-time is the product of the four-dimensional 
Minkowski space A^4 (the brane) and a d-dimensional 
torus T'^ and treats gravity in the linear approxi- 
mation gMN = Vmn + xohMN iVMN has mostly 
minuses). To calculate gravitational bremsstrahlung 
classically one has to extend the ADD setup beyond 
the linearized level. For this, one expands the metric 
further as Hmn — *■ hMN + xdSHmn and adds the cu- 
bic interaction terms to the Fierz-Pauli lagrangian. 
Equivalently, as in the standard theory of gravita- 
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tional radiation in four dimensions, one may expand 
the D-dimensional Einstein tensor up to quadratic 
terms: 

Gmn = -^□('/'MAT + xoSiiMN) - ^Xd^mn, (2) 

where the last term plays the role of a gravitational 
stress-tensor, whose form is dimension independent. 
Here h = ■q'^^hMN, ipMN = hMN - r]MNh/2 is the 
trace-reversed metric perturbation, □ = rj^^ dAjdN 
is the /^-dimensional d' Alembertian and the coupling 
constant is defined by = WttGu. Symmetriza- 
tion and alternation over indices is understood with 
1/2, while raising/lowering of indices is meant with 
the flat metric. 

Since Z^-dimensional Einstein equations imply the 
_D-dimensional Bianchi identities, which in turn im- 
ply the £)-dimensional geodesic equations for the par- 
ticles, going beyond the linearized theory might con- 
tradict the assumption of matter confinement on the 
brane. However, it turns out that it is enough to as- 
sume that to zeroth order in xjj particles move on the 
brane. Then, the corresponding zeroth order stress- 
tensor also lies on the brane 



(3) 



Its first order perturbation STmn, due to the lin- 
earized gravitational interaction, is also confined on 
the brane. Indeed, to first order the wave equation in 
the flat harmonic gauge Om^P^^^ = is 
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The source term is constructed neglecting gravity, 
hence it is flat-space divergenceless too. Since parti- 
cles move freely in this order, this equation describes 
non-radiative Lorentz-contracted gravitational poten- 
tials. Radiation appears in second order in >cd and is 
described by the field SipMN satisfying (again in the 
flat harmonic gauge) the equation 

aSipMN = —^dtain, tmn = STmn + Smn, (5) 

where S mn is quadratic in the first order gravitational 
potentials. The right hand side of ^ is divergenceless 
by virtue of the expanded Einstein equations. There- 
fore, the effective source of radiation is the sum of the 
perturbation STmn of the matter tensor, caused by 
the first order gravitational interaction and the gravi- 
tational stress tensor Smn, constructed from the first 
order metric perturbations. This tensor is not con- 
fined on the brane, but extends to the bulk. Note 
that STmn and Smn are not separately gauge invari- 
ant. Thus, beyond the linearized level the effective 
source of radiation is no more confined on the brane. 

Denote the I?-dimensional coordinates as x^'^ — 
(x^, y^), where a;^, /i = 0, 1, 2, 3, lie on the brane and 

, k = I, . . . ,d label the points of T''. Imposing peri- 
odicity conditions hMN{x,y''+2TTR) = hMN{x,y^) we 
obtain an infinite number of four-dimensional massive 
modes 

hMN{x'') = ^ /zX,^(x)e»'=^'"/« (6) 

with masses m^ — q^, where qk — rik/R is the 
quantized momentum on the torus, and which can 
be regrouped into spin 2,1,0 massive fields d, [l6j . 



But, being interested in the total radiation in all 
modes, it is more convenient to think about the 
metric perturbation as a D-dimensional massless 
field with discrete momenta in the extra dimensions. 
The four-dimensional fields h\ij^{x^) are further ex- 
panded into Fourier integrals defined by '^{x) = 
/^'(g)e-*«-^d4g/(27r)4. 

Consider the small-angle collision of two point 
masses m, m'. Assuming that to zeroth order in xd 



their world-lines lie on the brane 



'iT),x'- 



(t'), one verifies that, as expected, they remain 
on the brane even after the gravitational interaction 



/A" 



is switched on. Therefore, 
z'^(t) =fe'^ + ^r + (5z^, z'^(t) = -f fc'^, (7) 

771 to' 

where p'^ and p''^ are momentum parameters. Al- 
though the true initial momenta = m lim i^(T) 

T — > — OO 

and P"^ ~ to' lim z''^{t) differ from p^^, p'f^, they 

stiU satisfy s ^ {P + P'f = [p + p'f [i3- It is 
convenient to work in the rest frame of one of the par- 
ticles, m', choosing the coordinate axes on the brane 
so that^ = m'(l, 0,0,0), p'' = m7(l,0,0,7;), 7 = 
1/ y\ — v^ . Also, with no loss of generality one may 
set 6^ = (0, 6, 0, 0) and h'^^ — 0, so that b is the impact 
parameter. Finally, one may think of brane localized 
vectors as D-dimensional vectors with zero bulk com- 
ponents, e.g. p*^ = (p'^, 0, . . . , 0). 

In the linearized theory the superposition princi- 
ple implies that both the zeroth and the first order 
energy-momentum tensor is the sum of the contribu- 



tions of the two masses, i.e. Tj 



MN 



=Tmn+ T MN and 



STmn —STmn+ ST mn- The stress term Smn, on 
the other hand, represents a collective contribution, 
which cannot be attributed to any one of the par- 
ticles. Solving the wave equation ([5]), one finds that 
the radiation amplitude consists of three terms, which 
correspond to radiation from masses and from their 
gravitational stresses. This is similar to the structure 
of the Born amplitude in quantum theory , which 
in the case without extra dimensions and in the low 
frequency limit is known to coincide with the classi- 
cal result. However, like in elastic scattering [l3| in 
ADD {d ^ 0), the classical treatment is essentially 
non-perturbative in the quantum sense. 

The spectral-angular distribution of the emitted en- 
ergy in gravitational bremsstrahlung during the colli- 
sion is given by 



dujdQ 27r2y 



J2Y.\^MN{k)e^U 



MN\2 



(8) 



nel,'' pol 



where 



TMNik) =STMN{k)+ ST MN{k) + SMNik) (9) 

is the Fourier transform of the right hand side of (O . 
The first two terms have only the brane components 
M,N = 0,1,2,3, while the third is truly multidi- 
mensional. Introduce next the massive wave vector 
fc^ = {y/uj'^ + m^, wn), with m^ = k^, m = rn/R. 
The three-dimensional unit vector n will be parame- 
terized by the spherical angles 9, ip^ in the usual way. 
Alternatively, one can think of the radiation wave vec- 
tor as the D-dimensional null vector k^^ = (fc'', k'). 
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The polarization tensors of the emitted radiation are 
chosen according to the Z3-dimensional picture. They 
are D{D — 3)/2 symmetric transverse traceless tensors 
orthogonal to k^' and to each other. To construct 
them start with D — A unit space-like mutually or- 
thogonal vectors e^, a = 3, 4, . . . , I? — 2, orthogonal 
also to p*^, p'*^, k^^ and 6^^. Their contractions with 

m m 

5TMN{k), ST A/Af (fc) and SMN{k) will be zero. Take 
in addition the two vectors 

k ■ pp'^ — k ■ p'p^^ + k'^^mm'vj cos 6 



e2M 



mv"f k ■ p' sin 9 

eMNLPA3...AD-2P^P'^k^<^3^---(^D~2^ 

mvj k ■ p' sin 9 



(10) 



satisfying 62 ■ p = £2 ■ p' — ei ■ p' = 0. Using those, one 
builds two chiral graviton polarizations, which belong 
to the brane: 

(11) 



e'J = (ef ±zef)/\/2. 



and a third one extending to the bulk 



.MN 



= M V ef -{D- 4)e 



(Af N'j 



(12) 



where A/" is a normalization factor. The remaining 
polarization tensors contain at least one vector from 
the set {e^} and give zero contraction with tm n ■ 

The computation of the bremsstrahlung radiation 
proceeds as follows. One finds the retarded fields gen- 
erated by the unperturbed particle trajectories and 
substitutes to the particle geodesic equations to get 
the first order corrections to the particles' motion. 
These are used to build the perturbations of the par- 
ticle energy-momentum tensors, whose Fourier trans- 
forms are the first two terms in the radiation source 
([9]). The result is expressed in terms of Macdonald 
functions with argument Wn — [w^ + m^&^]^/^, wq = 
k ■ pb/lrwyv) and reads 



6Tf,,{k)=- 



E 



. Ki{W n) 

(wo) 



2 f^'^ 



+2 



j-KQ{'Wnn+ -^is:o(m„) 7— -1 



Wq 



, (13) 



where w'q — k ■ p' b/{'m'jv), cr^^ = [p^p^^k ■ b — 2k ■ 
pp(^6^-)]/m^ and K\{w) = K\{w). The second 

ni 

term 5T ^j^^{k) is obtained from by the substitu- 
tion {in,m' ,Pf^,p'^) — > {m' ,m,p'^,p^) and e^**^''^ 1. 

Next, one projects over polarizations. Only the first 
three contribute in the chosen gauge, and all of them 
have zero contractions with the m' term in Q . So, the 
total amplitude receives two contributions: (a) from 
the moving mass ([T3| , and (b) from the stress term in 
(l9|), whose significance will be discussed shortly. 

The next step is to sum over the interaction modes. 
Assuming that a large number of modes is excited, one 
can replace the summation over by integration: 



E 



/n 



Jfnd'n^Y^Jf,„qt'dq„, (14) 



where 5^ = n^/R. Use was made of the fact that 
the summand depends only on the absolute value of 
n (m^ = nV-R^)- The result is {l^: 



E ^^a(u;) 



{2nRyb^f K^+,/^iwo). (15) 
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Summation over modes in the amplitude is the clas- 
sical counterpart of integration over transverse mo- 
menta of virtual gravitons, which leads to tree-level 
divergences in ADD However, classically the re- 
sult is finite like in elastic scattering [iTj . 

Approximating summation over transverse modes 
with integration leads effectively to the situation in 
uncompactified Minkowski space Md- There too, the 
resulting particle source is expressed in terms of the 
Macdonald functions, now with argument wq. The 
properties of the Macdonald functions imply that the 
corresponding radiation amplitude is concentrated in 
a narrow cone 9 < I/7 with a frequency cutoff at 
ujc 27^/6. This is also what typically happens in 
electromagnetic bremsstrahlung in flat space. How- 
ever, the spectral-angular distribution in the case of 
gravitational interaction is very different, because of 
the special role of Smn- In 4d gravity it was shown 
that the contributions of (fT3|) and 5*^,^ cancel in the 
above spectral-angular region [l^. Physically, this is 
due to the fact that the radiation, emitted by particles 
following ultrarelativistic time- like geodesies, follows 
null geodesies, which are close to the former and give 
an effective formation length of radiation in a given 
direction 7 times larger than in flat space p^ . 

No surprise, the situation is similar in x T*^. 
The contribution from the stresses n , which is the 
classical counterpart of the amplitude involving the 
three-graviton vertex, is rather complicated. How- 
ever, it can be shown [l3| that like for d = the lead- 
ing term in the S'-amplitude for 7 3> 1 exactly can- 
cels (IT^ . The non- leading term, integrated over the 
transverse virtual momenta, is also expressed in terms 
of Macdonald functions, whose argument w'q ~ cub/j 
does not depend on angles. Therefore, the radiation 
does not exhibit sharp anisotropy and the frequency 
cutoff (from the condition Wq ~ 1) is w ^ uJct — j/b. 
Projection on e± leads to (£3 gives a similar contribu- 
tion): 



(cos 2ip q= 2i) cos 9' Kii/2 



4(27r)d/2+i 

{i cos^ Lp ^ sin 2(^)i^(;/2+i ~sin2(p(i sin(p ± \)Kdl2 



where ^ = 1 — w cos 9 and 9' is the radiation polar 
angle in the rest frame of m (sin0' = sin 0/7^). 

Finally, upon integration of ([5]) over O and w one 
obtains the radiated energy 



'2 ,^6 



mm K 



D^d+3 



(16) 



with a known dimension-dependent coefflcient. Qual- 
itatively the dependence on b and 7 can be under- 
stood as follows. The averaged over emission an- 
gles amplitude-squared is (|t±P) ~ 6^^*^, and i?rad ^ 
(kitP)^^ A/'cff, where AfeB ~ {Ri/bY is the effective 
number of emitted light modes [14l |. 

Relation (|16p was obtained in the rest frame of par- 
ticle m' . To pass to the CM frame, one calculates the 
relative energy loss (radiation efficiency) e = i?rad/£', 
and expresses the result in terms of the Lorentz factor 
in the CM frame via (for m — m') j'^^ = (1 + 7)/2: 
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/rc \ 3(<i+l) 

e = (17) 

The two new features of p7|) are: (a) the large factor 
7^^"'"^ due to the large number of light KK modes 
involved both in the gravitational force and in the 
radiation, and (b) the growing with d coefficient: 



d 


1 


2 


3 


4 


5 


6 




10.1 


184 


3359 


6 • 10^ 


1.06- 10^ 


1.8-10^ 


rs 


3.45 


1.88 


1.46 


1.29 


1.21 


1.17 


be 


196 


7.90 


3.15 


2.11 


1.72 


1.53 



The Table has rs and be in TeV ^ evaluated for 
~ ITeV and ^ ~ 14TeV. The classical descrip- 
tion of small angle ultrarelativistic scattering is valid 
for impact parameter in the region rs < b < be , where 
be = TT-i/^ [rid/2)GDs/hc^f ^ rsirs/Xuf'" is 
the scale beyond which (for d 7^ 0) the classical no- 
tion of particle trajectory is lost [6]. Another restric- 
tion comes from the quantum bound on the radiation 
frequency fuoct < rn"/, which is equivalent to 6 > Ac = 
h/ (mc). For d ^ the two conditions overlap provided 
Ac < be- To estimate e set & Ac in p7|) to obtain 
e = Bd{sm/MlY+^ [Bd = 7.4, .8, .6, .9, 1.9, 5.6, 21 for 
d = 0, 1, . . . , 6). Thus, a simple condition for strong 
radiation damping is 

sm > Ml , (18) 

which may well hold for heavy particles and nuclei 
with LHC energies and cosmic rays. For example, for 



-y/s = 14 TeV and m ~ .2 TeV all conditions are met 
for d — 1,2, and at the quantum boundary b = Ac 
one has ei ~ 5 x lO'', £2 — 10®. For protons in LHC, 
Ac > be, and (fTT]) does not apply for d^ 0. For d = 
and = Ac it gives e = .25. 

To conclude, our analysis shows that gravitational 
bremsstrahlung may lead to strong radiation damping 
in transplanckian collisions. Therefore, (a) one may 
have to include the reaction force [lB| in the study 
of BH production, which might even exclude the for- 
mation of a CTS. Incidentally, there are indications 
that gravitational collapse of an oscillating string does 
not take place, once gravitational radiation is taken 
into account [20|. (b) Our results also imply that 
bremsstrahlung is a strong process leading to missing 
energy signatures in transplanckian collisions, which 
may further constrain the ADD parameters [14| . 
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